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ABSTRACT 



We study the effect that non-equilibrium chemistry in dynamical models of 
collapsing molecular cloud cores has on measurements of the magnetic field in 
these cores, the degree of ionization, and the mean molecular weight of ions. 
We find that OH and CN, usually used in Zeeman observations of the line-of- 
sight magnetic field, have an abundance that decreases toward the center of the 
core much faster than the density increases. As a result, Zeeman observations 
tend to sample the outer layers of the core and consistently underestimate the 
core magnetic field. The degree of ionization follows a complicated dependence 
on the number density at central densities up to 10 5 cm -3 for magnetic models 
and 10 6 cm -3 in non-magnetic models. At higher central densities the scaling 
approaches a power-law with a slope of -0.6 and a normalization which depends 
on the cosmic-ray ionization rate £ and the temperature T as (£T) 1//2 . The mean 
molecular weight of ions is systematically lower than the usually assumed value 
of 20 — 30, and, at high densities, approaches a value of 3 due to the asymptotic 
dominance of the Hjj~ ion. This significantly lower value implies that ambipolar 
diffusion operates faster. 

Subject headings: ISM: molecules - ISM: clouds - ISM: dust - magnetic fields - 
MHD - stars: formation - ISM: abundances 



1. Introduction 

Prestellar cores are the earliest identified phases of star formation. They are centrally 
condensed, likely to be gravitationally bound and destined to ultimately form stars or clus- 
ters. They have power- law density pro files with flat inner regions, implying that no proto - 



star has yet formed in their center (e.g. IWard-Thompson et al.lll999t iBacmann et al.ll2000l ). 



1 Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA 
2 Current address: Max-Planck Institut fur Radioastronomie, 53121 Bonn, Germany 



-2 - 



Their linewidths are thermalized, implying no significant turbulence sup port (e.g. iMyers 



1983; Barranco fc Goodman 1998; Goodman et al 



1998 



Kirk et al.l 120071). They show no 



evidence for substantial a mounts of rotation (IGoodman et al.l Il993l ; iBarranco &: Goodman 
19981 ; ICaselli et al.l l2002al ) . Intensity r atios of different transitions sho w that starless cores 



are almost isothermal, with T ~ 10 K (Tafalla et al 



1998 



2002 



20041 ). although there may 



be small variations of the order of a few K (e.g. lEvans et al.ll200ll ). Their evolution depends 
critically on the physical conditions of the star forming region in which they are embedded. 
In addition, their relatively quiescent dynamical state allows for an easier interpretation of 
observations because of the absence of the effects of thermal and dynamical feedback from a 
central protostar. For these reasons prestellar cores are ideal probes of the ultimate stellar 
origins: the initial conditions of star formation. 



Prestellar cores are observed to be magnetic (see, e.g. iHeiles fc Crutcherl 120051 . for a 
review). The magnetic field can provide support against the self-gravity of a molecular 
cloud or core, and thus affect its dynamical evolution. The amount of magnetic support 
is quantified by the mass to magnetic flux ratio M/$b of the object under consideration. 
There is a critical value for the mass-to-flux ratio, 



63G 



1/2 



(i; 



(IMouschovias &: Spitzerl Il976l ) where G is the gravitational constant. If the mass-to-flux 
ratio of an object exceeds this value (the object is magnetically supercritical), then the 
magnetic field is not strong enough to support the object against its own self-gravity, and 
the object contracts dynamically. The opposite is true for mass-to-flux ratios below the 
critical value (magnetically subcritical objects). Magnetic models of cloud fragmentation 
and core collapse predict that prestellar cores are magnetically supercritical, dynamically 
collapsing fragments, formed in magnetic parent clouds. These parent clouds which can be 
magnetically subcritical as a whole, and the supercritical fragments are formed through the 
process of ambipolar diffusion. The latter is the process of neutral particle diffusion through 
ions and magnetic fie ld lines toward centers of grav ity which increases the mass-to-flux ratio 



of the fragment (e.g. iFiedler &: Mouschoviaslll993l ). 



Determining the amount of magnetic support in a molecular cloud as a whole and in 
individual molecular cloud cores is made complicated by at least the following three effects. 

First of all, Zeeman measurements only trace the component of the magnetic field along 
the line of sight, and for this reason geometrical considerations due to unknown cloud and 
core orient ations limit the amount of information that can be obtained on an object -to-object 
basis (e.g. IShu et al.l 1 19991 ; iTroland fc Crutcherl 120081 ; IMouschovias fc Tassisl 120091 ). 
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Second, in the case of Zeeman measurements, the magnetic field strength is convolved 
with the abundance profile of the tracer molecule, and, as a result, it is non-trivial to de- 
termine which parts of a core contribute most to a finite-beam Zeeman measurement of its 
magnetic field. Understanding the time evolution and spatial variation of the abundance of 
different species is essential. Observations show that when such measurements are made us- 
ing different molecular lines, they reveal different values of the magnetic field which cannot 
be reconciled. A characteristic exa mple is CN and OH Z eeman obse rvations of the same 
object s yielding different results (e.g. iFalgarone et al.ll2008l ). Similarly. ICrutcher fc Troland 
( I2000l ) have argued, in the case of L1544, that OH data (and the associated Zeeman measure- 
m ent of the magn e tic fie ld) do not sample the small, dense core observed in dust emission, 



so ICrutcher et al.l (120041 ) concluded that their Chandrasekhar- Fermi-measured plane-of-sky 
magnetic field may be disc repant from the Zeeman-measured line-of-sight magnetic field of 
Crutcher &: Trolandl ( 120001 ) . This may be a result not only of geometrical projection effects, 
but also because different measurements correspond to different parts of the core. To verify 
the origin of such discrepancies, we need to understand how different molecules trace different 
parts of cores and thus different magnetic field values. This becomes especially important 
when attempting to study the variation of the magnetic field strength with density (e.g. 
Crutcher et aDboiol ). 



Third, by the time the column density contrast between core and cloud becomes large 
enough for the core to be detected with high statistical significance and to be studied in 
detail, the core has already left the quasistatic contraction phase and has entered the phase 
of dynamical collapse, even in models which are originally heavily magnetically supported 
( jTassis fc Mouschoviasl 12004 ). At this stage, the mass-to-flux ratio is alrea dy larger than 
critical regardless of its value at larger scales (e.g. iMouschovias et al.ll2006l and references 
therein). At even higher central densities, during the very advanced stages of collapse, the 
dynamical evolution of cores and the resulting density and velocity pr ofiles have been shown 
to be quite insensitive to the initial value of the mass-to-flux ratio ( jTassis fc Mouschovias 
2007bh . 



In addition to the mass-to-flux ratio, the effect of the magnetic field on the dynamics of 
molecular clouds and molecular cloud cores also depends on the degree of ionization of the 
cloud or core. The ambipolar diffusion timescale is proportional to the degree of ionization, 

71; 



^ad oc 



72 H o 



(2) 



( jCiolek fc Mouschoviaslll993[ ). This dependence is physically straight-forward to understand: 
a higher degree of ionization implies a better coupling between magnetic field and matter, 
which results in an increased resistance of the ions as the neutrals drift past them and con- 
sequently an increased ambipolar diffusion timescale. The degree of ionization of a cloud or 
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core is determined by the ionization rate (in dense cores the dominant ionization mechanism 
is cosmic ray ionization, and hence the relevant quantity is the cosmic ray ionization rate () 
and by the relevant recombination reactions. Because of its feedback effect on the dynamics, 
and especially in the case where the chemistry is out of equilibrium, the degree of ionization 
is non-trivial to obtain; ideally, it has to be self-consistently calculated through a model 
following both dynamics and chemistry simultaneously. 

For these reasons, the role of magnetic fields in the fragmentation of molecular clouds 
and the core formation and evolution process remains observationally uncertain and a hotly 
debated subject in the field. 



In a companion paper (ITassis et al.ll201ll . hereafter Paper I) we discussed our models of 



evolving prestellar cores which couple non-equlibrium chemistry with dynamics. Our exten- 
sive parameter study included magnetic and non-magnetic dynamical models with varying 
initial values of the mass-to- flux ratio or collapse retardation times respectively, varying C/O 
ratios, cosmic-ray ionization rates, and temperatures. In Paper I we focused on the evolu- 
tion of the abundances of the most abundant and commonly observed molecules, and their 
dependence on the various model parameters. 

Here, we focus on the properties of the model cores that can help elucidate the role of 
magnetic fields in core formation and evolution: the degree of ionization, the most common 
ions and the mean ion molecular weight, abundance profiles of molecules used in Zeeman 
observations, and how the magnetic field that would be measured through Zeeman obser- 
vations in our model cores depends on the actual strength of the magnetic field at the core 
center. 

This paper is organized as follows. A brief review of the different models we consider 
is given in §2J The results and their dependence on the various parameters we have studied 
are presented in ^J3j We summarize and discuss our conclusions in §H 



2. Models 

Since in this paper we are interested not only in the measurement of the magnetic fields 
but also in the evolution of the degree of ionization, we consider both magnetic and non- 
magnetic models. The details of the dynamical and chemical models of each class is discussed 
in detail in §2 and 3 of Paper I. Here, we very briefly review the parameters we have varied 
for each class of models presented and discuss in the following sections. These parameters 
include the temperature, the C/O ratio, the cosmic ray ionization rate, and a parameter 
controlling the time available for chemical evolution. The latter is the mass-to-flux ratio 
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Magnetic runs Non-magnetic runs 




Fig. 1. — Line types and colors used to denote each of the models studied in this work, unless 
otherwise noted. Solid normal-thickness red line: "reference" magnetic model; solid normal- 
thickness blue line: "reference" non-magnetic model. Dotted lines: "fast" models; dashed 
lines: "slow" models. Brown/purple lines: magnetic/nonmagnetic models with tempera- 
tures differing from the "reference" models. Orange/cyan shaded areas: variation in C/O 
ratio. Thin/thick solid red/blue lines: lower/higher cosmic-ray ionization rate magnetic/non- 
magnetic models (see text for details). 
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in the case of magnetic models, and the collapse delay time in the case of non-magnetic 
models (an initial time period during which chemistry evolves but the core does not evolve 
dynamically, representing an early stage of support due to turbulence which later decays). 

Our "reference" magnetic model has a mass-to-flux ratio equal to the critical value 
for collapse, a temperature of 10 K, a C/O ratio of 0.4, and a cosmic ray ionization rate of 
( = 1.3 x 10~ 17 s _1 . Our "reference" non-magnetic model has a collapse delay time (hereafter 
"delay") of lMyr, and values for the C/O ratio, temperature, and ( identical to those of the 
"reference" magnetic model. 

For magnetic models we examine two additional values of the initial mass to magnetic 
flux ratio: 1.3 times the critical value (a faster-evolving, magnetically supercritical model), 
and 0.7 of the critical value (a slower, magnetically subcritical model). For non-magnetic 
models we examine two additional values of delay: zero, and lOMyr. 

For each of these six dynamical models the carbon-to-oxygen ratio is varied from its 
"reference" value by keeping the abundance of C constant and changing that of O. The two 
other values of C/O ratio examined are 1 and 1.2. We have studied in this way a total of 18 
different models (9 magnetic and 9 non-magnetic). 

In addition, to test the effect of the temperature, we have varied each of the six basic 
dynamical models by changing T by a factor of ~ 1.5 from its reference value of 10 K and 
examined models with T = 7 K and T = 15 K. We have thus studied 12 models (6 magnetic 
and 6 non-magnetic) with temperature varied from its reference value. 

Finally, to test the effect of the cosmic ray ionization rate, we have studied four ad- 
ditional models (two magnetic and two non-magnetic), which have a "reference" value for 
the temperature, C/O ratio, and mass-to-flux ratio or delay (for magnetic and non- magnetic 
models respectively), but for which ( is varied by a factor of four above (£ = 5.2 x 10~ 17 s ) 
and below (( = 3.3 x 10~ 18 s -1 ) i ts "reference" value (covering the range of observational 



estimates (e.g. McCall et al.ll2003l : iHezareh et al.ll2008f ) 



These additional models bring the total of different models we have run and examined 
to 32. Figured] visually depicts these different models and the line type/color used to denote 
each one (unless explicitly noted otherwise). 
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Fig. 2. — Left column: upper panel: radial profiles of the z-component (solid lines) and 
r-component (dotted lines) of the B — field at four different times t\ — £4 (corresponding 
to the same snapshots for which the average B z traced by OH and CN is shown in Fig. 0}; 
lower panel: radial profiles of the number density at the same times t\ — t±. Right column: 
radial profiles of the OH (upper panel) and CN (lower panel) abundance taken at the same 
times as in the left column. 
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Fig. 3. — Column density of CO, NH 3 , OH, and CN plotted against the fractional radius 
of the core, for a central density of 10 6 cm -3 (corresponding to snapshot t% in Fig. [2]). The 
three red lines correspond to the different magnetic dynamical models (magnetically critical, 
subcritical and supercritical for the solid, dashed, and dotted lines respectively) for a refer- 
ence value of the C/O ratio equal to 0.4, while the yellow line corresponds to a magnetically 
critical model with C/O ratio equal to 1 (see Fig. [T|). 
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Fig. 4. — Mean value of the z-component of the magnetic field of a core traced by OH (solid 
lines) and CN (dotted lines) for the initially magnetically subcritical (diamonds), critical 
(x), and supercritical (+) cloud models. For comparison, the initial magnetic field of the 
model cloud is B z =4.2, B z =5.6, and B z =7.5 for the magnetically supercritical, critical 
and subcritical cloud models, respectively. 
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Fig. 5. — Evolution of central electron abundance versus central number density in magnetic 
and non-magnetic models. In all panels, red lines correspond to magnetic models with refer- 
ence values for C/O and temperature (dashed: magnetically subcritical; solid: magnetically 
critical; dotted: magnetically supercritical), and blue lines to non- magnetic models with 
reference values for C/O and temperature (dashed: lOMyr delay; solid: lMyr delay; dotted: 
no delay). The black dashed line corresponds to Eq. (J3J). Left panel: effect of a varying C/O 
ratio. The cyan and orange shaded areas (for non-magnetic and magnetic models respec- 
tively) correspond to a range of C/O values between 0.4 (reference value) and 1.2. Middle 
panel: effect of a varying core temperature. The thin and thick brown (purple) solid lines 
correspond to temperature of 7 K and 15 K respectively for magnetic (non-magnetic) models; 
the "reference" value for the temperature is 10 K. Right panel: effect of a varying cosmic-ray 
ionization rate. The thick and thin red (blue) solid lines correspond to ( a factor of four 
above and below the "reference" value for magnetic (non- magnetic) models (see Fig. [Q. 



- 11 - 



3. Results 

3.1. Zeeman- Traced Core Magnetic Field 

In this section we examine how the magnetic field value that is measured through 
finite-beam Zeeman observations compares to the magnetic field at the center of a core. 
We consider two molecules commonly used for Zeeman observations, OH and CN. Figure [2] 
shows radial profiles of the magnetic field (left column, upper panel; z— and r— components 
are shown with solid and dotted lines respectively), the OH abundance (right column, upper 
panel) and the CN abundance (right column, lower panel) at four time instances t\ — t 4 , 
each of which represents an increase in the H2 number density at the center of the core 
of an order of magnitude. The central number density profiles at these times are shown 
in the lower panel of the left column. These radial profiles correspond to the magnetically 
critical ("reference") dynamical model. Note that although CN and OH appear significantly 
de pleted at the very centr al parts of the core, these results are consistent with observations 



by iHily-Blant et al.l (120101 ) . who find that significant amounts of CN remain at densities of 
~ 3 x 10 4 cm -3 , where CO has already depleted. Comparing with the lower-left panel of 
Fig. El we find that significant depletion of CN (abundance more than an order of magnitude 
smaller than in the outskirts of the core) only sets in at densities higher than ~ 10 5 cm -3 . 
To facilitate comparison with such observations, we show, in Fig. [31 the column density 
of four molecules of interest (CO, NH 3 , OH, and CN) as a function of (linearly plotted) 
fractional radius of the core, for a central density of 10 6 cm -3 , corresponding to snapshot 
£3 in Fig. [21 when the central abundance of CN for example is more than five orders of 
magnitude lower than that in the outskirts of the core. The three red lines correspond to 
the different magnetic dynamical models (magnetically critical, subcritical and supercritical 
for the solid, dashed, and dotted lines respectively) for C/O ratio equal to 0.4, while the 
yellow line corresponds to a magnetically critical model with C/O ratio equal to 1. A peak 
at relatively large radius and then a plateau/mild decrease of the column density toward 
the center is a sign of depletion. In this context, CO and NH3 represent extreme cases of 
significant and little depletion respectively. OH and CN are intermediate cases, bot h peaked 



at smaller radii than CO, in agreement with the findings of IHily-Blant et al.l (120101 ) . If C/O 
is 1 rather than the fiducial 0.4 then CN exhibits a central decrease/plateau at much smaller 
radii, while CO and OH are significantly depleted. 

In Fig. H] the mean magnetic field that would be seen through Zeeman measurements 
by a O.lpc beam (extending out to a radius where the neutral number density falls to about 
10 4 cm -3 ) is plotted against the actual central value of the z— component of the magnetic 
field. The mock Zeeman-obtained "observed" value of the field is derived as follows. We 
calculate the average ^-component magnetic field in the magnetic core under consideration, 
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weighted by the number density of the corresponding molecule (OH or CN respectively) 
within a radial extent of 0.1 pc (the assumed size of the beam) at four different time instances. 
The core is assumed to be viewed face-on (down the z — axis). In this way, geometrical 
effects do not enter our calculation of the mock Zeeman observation (the z— component of 
the magnetic field, which we examine here, is oriented exactly along the line of sight), and 
we can instead focus only on the effect of depletion. 

We then compare each of the averages derived in this way to the value of the z-component 
of the magnetic field at the center of the core at the same time. In practice, this means that 
(B z ) is calculated through weighting by the OH or the CN number density in each radial an- 
nulus within the beam, respectively, assuming optically thin lines. Three different dynamical 
models are shown: magnetically subcritical (diamonds), critical (x) and supercritical (+). 
Solid lines are used for the magnetic field traced by OH and dashed lines for the magnetic 
field traced by CN. 

In all cases we can see that, although at low central densities the value of the Zeeman- 
traced magnetic field is close to the actual value of the central B z , at higher densities the 
measured field increases only mildly and shows an overall increase by only a factor of two, 
whereas the actual central value of B z has increased by an order of magnitude. Most of the 
increase in the Zeeman-measured field strength takes place at low densities, while at higher 
densities the measured field tends to saturate. The reason for this behavior is seen in Fig. |2j 
the abundance of OH and CN is falling at the center of the core much more rapidly with 
time than the density increases. As a result, the Zeeman-traced magnetic field corresponds 
to the outer layers of the core "onion skin" , where OH and CN have still relatively higher 
abundances 



3.2. Degree of ionization 

Figure shows the evolution of electron abundance (degree of ionization) versus central 
number density in magnetic and non-magnetic models. In all panels red lines correspond 
to magnetic models with reference values for C/O and temperature (dashed: magnetically 
subcritical; solid: magnetically critical; dotted: magnetically supercritical), and blue lines to 
non-magnetic models with reference values for C/O and temperature (dashed: lOMyr delay; 
solid: lMyr delay; dotted: no delay; see also Fig. [TJ). The left panel shows the effect of 
varying the C/O ratio, the middle panel the effect of varying the temperature, and the right 
panel the effect of varying the degree of ionization. For comparison we overplot with the 
black dashed line the usually adopted scaling [obtained through fits to earlier calculations, 
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Fig. 6. — Degree of ionization as a function of central number density n c for the three 
magnetic (upper panels) and non-magnetic (lower panels) models, with "reference" values 
for the C/O ratio and the ionization rate (. Left column: "fast" models; middle column: 
"reference" models; right column: "slow" models. Pie charts show the relative contribution 
of the dominant ions at the central densities denoted by stars. 
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e.g., lElmegreenl ( 119791 ); iNakand (119791 )] between degree of ionization and number density. 

Hi „ / n \ fe-i 
n 



10 5 cm -3 



(3) 



(IBasu fc Mouschoviasl 1 1 9941 ). with k = 0.5 and K = 3 x 10" 



It is obvious from Fig. [5] that at low central densities the out-of-equilibrium chemistry 
causes a variety of transitional effects, and the scaling of the degree of ionization with 
density cannot be expressed in a simple form such as Eq. (jHJ). For central densities higher 
than ~ 10 5 cm -3 for magnetic models and ~ 10 6 cm -3 for non-magnetic models the scaling 
does asymptotically approach a simple power-law. The magnetic models approach power-law 
scaling at lower central densities because of the increased time available to the chemistry to 
overcome transitional effects before a specific density is reached. The slope of this scaling 
(-0.6, corresponding to k — 0.4) is steeper than t he va lue usually adopted, but well within 
the uncertainties quoted in lBasu fc Mouschoviasl (119941 ). who give a ra nge for k between 0.3 
and 0.5. Our value for the slope is also consistent with the findings of ICaselh et all (j2002bl ) 
for the case of L1544. The normalization of our models is also consistently lower than Eq. 
([2]), but the discrepancy is not larger than the uncertainty on K (an order of magnitude). 

Except for the earliest evolutionary times that are dominated by transitional effects 
in chemistry, the slope of the scaling initi ally steepens (the value of k i n Eq. [3] decreases) 
with increasing density, in agreement with ICiolek &: Mouschoviasl ( 119941 ) . However, once k 
reaches 0.4, the scaling of the degree of ionization with density becomes a power law, in 
contrast to ICiolek fc Mouschoviasl ( 119941 ). who found that k continues to decrease and the 



scaling cont inues to steepen. This is like ly a result of the (equilibrium) chemistry network 
adopted by ICiolek fc Mouschoviasl ( 119941 ) , which was considerably more simplified than the 
(non-equilibrium) chemical model we use here, and which, as we will see below, resulted in 
a considerably different ion population, dominated by different species. 

The non-magnetic models show a qualitatively different behavior, with the scaling of 
the degree of ionization with density becoming more shallow before it later steepens. At 
sufficiently high central densities the non-magnetic models asymptotically approach the same 
power laws as the magnetic models with the same value of £. The normalization of the scaling 
changes with C approxim ately as £ 1//2 , consistent with the usually assumed scaling (see, e.g., 
McKee fc Ostrikerl 120071 Eq. 24). The dependence at high densities on temperature is small 
(the asymptotic normalization being oc T 1//2 ), however at low densities it can be considerable. 
Similarly, the dependence on the C/O ratio can be appreciable at low densities. 

At sufficiently late evolutionary times and their associated high central densities mag- 
netic and non-magnetic models with the same value of ( and T converge to the same scaling, 
regardless of delay time, C/O ratio, or mass-to-flux ratio. We note that the high-temperature 
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non-magnetic model has not yet converged to its asymptotic form by the highest density dis- 
played in Fig. |5) for this reason, we have extended this run to even higher central densities, 
and we have confirmed that above 10 7 cm -3 the degree of ionization does indeed approach 
the corresponding scaling of the magnetic model with the same values of ( and T. 



3.3. Dominant ions and mean molecular weight 



Figure [6] shows the degree of ionization at the center of the core as a function of central 
number density n for the three magnetic (upper panels) and non-magnetic (lower panels) 
models, with "reference" values for the C/O ratio, temperature, and cosmic-ray ionization 
rate (. "Reference" dynamical models are shown in the middle column, with "fast" and 
"slow" models in the left and right columns, respectively. The pie charts above the lines 
show the relative contribution of the dominant ions to the total ion population at the densities 
marked by stars. 

Due to our adopted initial condition, which assumes all C to be ionized, C + is initially 
the dominant ion. However, once the chemistry starts operating, this evolves to a different 
mixture. Non-magnetic models have a more diverse mixture of ions at early times, with 
significant contributions from HCO + and H 3 + . 

At late evolutionary times Hjj~ is the dominant ion in all models, with H + and He + being 
the most important secondary contributors (grain gro wth, not considered here, may alter the 
relative importance of Hjj~ and H + in fav or of H + . seelFlower etal.l 120051). This is an impor 



t ant difference from previous stud i es fe.g.lCiolek fc Mouschoviaslll993l ; iDesch fc Mouschovias 



200ll ; iTassis fc Mouschoviasll2007l ; iMcKee et al.ll2010l ). which have assumed that the ion pop- 
ulation at comparable densities is dominated by HCO + or by metallic ions such as Na + , both 




Fig. 7. — Mean molecular weight /ij at the center of the core, plotted as a function of central 
number density n c . Panels, lines, and colors as in Fig. [5] (see also Fig. [I]). 
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of which have much higher molecular weight, implying that the ion fluid has much higher 
inertia at a given number density. Interestingly, the highest-density pie chart for all magnetic 
models and the "slow" non-magnetic model are very similar, so this ion distribution appears 
to be the one to which the ionization profile of the core settles at late enough times. 

For magnetic models the ion with the next highest contribution at low densities is C + . 
Non-magnetic models, on the other hand, feature additional substantially contributing ions, 
such as C 3 H+ H 3 0+, HCO+. 

Figure [7] shows the evolution with central number density of the mean molecular weight 
of ions /ij. Panels, lines and colors correspond to various magnetic and non-magnetic models 
as described in Fig. [5] (see also Fig. [TJ). 

As expected from the results of Fig. contrary to the usually adopted values of ^ which 
are in the range of 20-30, the mean molecular weight of ions generally stays below 20. At 
high densities \ii asymptotically approaches 3 for all cases except the high-tem perature non 



magne tic model, consistent with Hjj" dominating the ions in these cases (see also lCaselli et al. 



2002b|). The non- magnetic "slow" model appears significantly different (has a much lower 
mean molecular weight of ions at low densities) than all other models because before the gas 
starts to collapse, the chemistry has already had lOMyr to evolve, and the mixture of ions 
is already close to the characteristic "late time" mix of the other models. 

The low values of (for most of the central number density range) in the magnetic 
models implies that the ambipolar diffusion timescale is shorter than the one computed using 
the usually adopted higher value for /ij, due to the smaller inertia of the ion fluid (which is 
attached to the magnetic field lines at these densities) and its associated decreased ability 
to provide resistance to the diffusion of the neutral fluid towards centers of gravity. 

The C/O ratio is shown to affect the mean molecular weight appreciably. Its effect is 
stronger in the non-magnetic models, since oxygen and carbon-bearing molecules are more 
abundant in the gas phase in faster-evolving models such as the "reference" and "fast" non- 
magnetic models. Indeed, the "slow" nonmagnetic model shows a much smaller sensitivity 
to the C/O ratio, comparable to that of the magnetic models. Similarly, the temperature has 
a significant effect on /Zj, mainly through its effect on the chemistry of N-bearing molecules 
(see also discussion in Paper I). Specifically, we have verified that the increase in the mean 
molecular weight with temperature in the non-magnetic model (which persists even at high 
densities) is due to an increased abundance of N 2 H + . By contrast in the case of the magnetic 
high-temperature models, H + is the dominant ion at intermediate central densities and \ii 
decreases below three. 
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4. Conclusions 

We have examined the effect that non-equilibrium chemistry in dynamical models of 
collapsing cores has on molecular abundances and thus on measurements of the magnetic 
field in these cores, on the degree of ionization and on the mean molecular weight of ions. 
We have considered both magnetic and non- magnetic models and models with different C/O 
ratios, cosmic-ray ionization rates, and core temperatures. 

We have found that molecules usually used in Zeeman observations of the line-of-sight 
magnetic field (OH and CN) have an abundance that decreases toward the center of the core 
much faster than the density increases (see Fig. [2]). Thus, Zeeman measurements tend to 
sample the outer layers of the core and consistently underestimate the core magnetic field, 
especially for higher- density cores with higher magnetic field strengths. 

The degree of ionization was found to follow a complicated dependence on the num- 
ber density in early evolutionary phases, when central densities are less than 10 5 cm -3 for 
magnetic models and 10 6 cm -3 for non-magnetic models. At higher densities the scaling 
approaches a power-law with a slope of -0.6 (slightly steeper than the usually assumed value 
of -0.5) and with a normalization that scales with the product of ( and T to the 1/2 power; 
however, we note that for non-magnetic models with temperatures higher than 10 K, the 
dependency does not attain its asymptotic form until much higher central densities. 

The mean molecular weight of the ions was found to be systematically lower than the 
usually assumed value of 20 — 30, and, at high densities, to asymptotically approach a value 
of 3 for all models, due to the asymptotic dominance of IPt. The only exceptions are the 
15 K non-magnetic models (no delay and lMyr delay), for which N 2 H + becomes dominant 
at higher densities due to the sensitivity of nitrogen chemistry to temperature (see paper 
I, §5.3). At low densities the dominating ions and the associated value of //j evolve with 
density, and exhibit sensitivity to the value of the C/O ratio, (, and the temperature. The 
considerably lower value of /z, for the magnetic models compared to the usually assumed one 
implies that ambipolar diffusion operates faster. 

We thank Paul Goldsmith and the anonymous referee for insightful and constructive 
comments that improved this paper. This work was carried out at the Jet Propulsion Lab- 
oratory, California Institute of Technology, under a contract with the National Aeronautics 
and Space Administration. ©2012. All rights reserved. 
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